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A Symbolic Summation Approach to Find 
Optimal Nested Sum Representations 

Carsten Schneider 



Abstract. We consider the following problem: Given a nested sum expres- 
sion, find a sum representation such that the nested depth is minimal. We 
obtain a symbolic summation framework that solves this problem for sums 
defined, e.g., over hypergeometric, g-hypergeometric or mixed hypergeometric 
expressions. Recently, our methods have found applications in quantum field 
theory. 



1. Introduction 

Karr's algorithm (KarSl; Kar85) based on his difference field theory provides 
a general framework for symbolic summation. E.g., his algorithm, or a simplified 
version presented in (Sch05c), covers summation over hypergeometric terms (Gos78; 
Zci9f), g-hypergeometric terms (PR97) or mixed hypergeometric terms (BP99). 
More generally, indefinite nested product-sum expressions can be represented in his 
nS-difference fields which cover as special cases, e.g., harmonic sums (BK99; Vcr99) 
or generalized nested harmonic sums (MUW02). 

In this article much emphasize is put on the problem how these indefinite nested 
product-sum expressions can be simplified in a ni]*-field. E.g., with our algorithms 
we shall compute for the sum expression 



(1.1) ^ = E — 

r=l 

the alternative representation 



1=1 



(1.2) B ^ ^(^H^, + 2H^ + 6iH^ + l)Hi'^H^ + 3{Hi^^f + i8Hn + 4)Hl'H6Hi^^) 

where A{n) = Bin) for all n > and where the nested depth of B is minimal; 

Hn = Y^'k=i i denotes the nth harmonic numbers and Hn'^ = X]fc=i p" are the 
generalized harmonic numbers of order o > 1. 
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In order to accomplish this task, we exploit a new difference field theory for 
symbolic summation (Sch08) that refines Karr's HS-fields to the so-called depth- 
optimal nE*-fields. In particular, we construct explicitly a difference ring monomor- 
phism (Sch09) which links elements from such a difference field to elements in the 
ring of sequences. Using this algorithmic machinery, we will derive for a given nested 
product-sum expression A a nested product-sum expression B with the following 
property: There is an explicit A G N = {0, 1, 2, . . . } such that 



and among all such alternative representation for A the depth of B is minimal. 

From an applicational point of view our algorithms are able to produce d'Alem- 
bertian solutions (N6r24; AP94; SchOl), a subclass of Liouvillian solutions (IIS99), 
of a given recurrence with minimal nested depth; for applications arising from 
particle physics see, e.g., (BBKS07; MS07; BBKS08; BKKS09a; BKKS09b). The 
presented algorithms are implemented in the summation package Sigma (Sch07), 
that can be executed in the computer algebra system Mathematica. 

The general structure of this article is as follows. In Section 2 we introduce 
the problem to find optimal sum representations which we supplement by concrete 
examples. In Section 3 we define depth-optimal ni]*-extensions and show how 
indefinite summation can be handled accordingly in such fields. After showing how 
generalized d'Alembertian extensions can be embedded in the ring of sequences in 
Section 4, we are ready to prove that our algorithms produce sum representations 
with optimal nested depth in Section 5. Applications arc presented in Section 6. 

2. The problem description for indefinite nested sum expressions 

Inspired by (BL82; NP97) one can consider the following general simplification 
problem. Let X be a set of expressions (i.e., terms of certain types), let K be a field"'^, 
and let ev : X X N — > K with {x,n) i— > x{n) be a function. Here one considers the so 
called evaluation function cv as a procedure that computes x{n) for a given x X 
and n G N in finite steps. In addition, we suppose that we are given a function 
t) : X — > N which measures the simplicity of the expressions in X] subsequently, we 
call such a triple (X,eY,d) also a (measured) sequence domain; cf. (NP97). 
In this setting the following problem can be stated: Given A ^ X; find B E X and 
A G N s.t. (1.3) and s.t. among all such possible solutions 5(i?) is minimal. 

In this article, the expressions X are given in terms of indefinite nested sums 
and products and the measurement of simplicity is given by the nested depth of 
the occurring sum- and product-quantifiers. Subsequently, we shall make this more 
precise. Let {X, ev, £)) be a measured sequence domain with ev : X x N ^ K and 
d:X^N. 

Example 2.1. Let X = K{x) be a rational function field and define for / = 2 ^ 
]K(a;) with p,q E K.[x] where q =^ and p, q being coprime the evaluation fimction 



Subsequently, all fields and rings are commutative and contain the rational numbers Q. 



(1.3) 



A{k) = B{k) Vfc > A 



(2.1) 




OPTIMAL NESTED SUM REPRESENTATIONS 



3 



here p(k),q{k) is the usual evaluation of polynomials at fc G N. In particular, we 
define = 1 if / e K(a:) \K and c)(/) = if / £ K. In the following {K{x), ev, O) 
is called the rational sequence domain. 

Example 2.2. Suppose that K = IfC'(qi, . . . , 5,„) is a rational function field 
extension over K' and consider the rational function field X := K(.t,xi, . . . 
over K. Then for f = ^ X with p,q (z K.[x, xi, . . . , Xm] where g ^ and p, q 
being coprimc we define 

ifq{k,q'[,...,q^^) = 
'iiq{k,q'[,.-.,qt)¥'0. 



(2.2) ev(/,A;) = 




Note that there is a 5 G N s.t. for all fc e N with k > 5 wc have q{k, q^, . . . , q^) ^ 0; 
for an algorithm that determines i5 see (BP99, Sec. 3.7). In particular, wc define 
D(/) = if / e K, and 5(/) = 1 if / ^ K. In the following (IK(a;, xi, . . . , Xm), ev, 5) 
is called q-mixed sequence domain. Note: If m ~ 0, we are back to the rational 
sequence domain. If we restrict to the setting K(a;i, . . . , Xm) which is free of cc, it 
is called q-rational sequence domain. 

More generally, X can contain hypergcomctric, q-hypergeometric or mixed hy- 
pergeometric terms; for instance see Example 4.10. Over such a set X wc con- 
sider the set of (indefinite nested) product-sum expressions denoted by ProdSum(X) 
which is defined as follows. Let © Sum ,Prod be operations with the signatures 

© : ProdSum(X) x ProdSum(X) ProdSum(X) 

® : ProdSum(X) x ProdSum(X) ProdSum(X) 

Sum : N x ProdSum(X) ProdSum(X) 

Prod : N X ProdSum(X) ProdSum(X). 

Then ProdSum(X) 3 X is the smallest set that satisfies the following rules: 

(1) For any f,g e ProdSum(X), f ®g € ProdSum(X) and f ®g & ProdSum(X). 

(2) For any / G ProdSum(X) and any r G N, Sum(r, /) G ProdSum(X) and 
Prod(r,/) G ProdSum(X). 

The set of all expressions in ProdSum(X) which are free of Prod is denoted by 
Sum(X). Sum(X) is called the set of (indefinite nested) sum expressions over X. 

Example 2.3. Given (X, ev,c)) from Example 2.1 with IK = Q and X = Q(a;) 
the following indefinite nested sum expressions are in Sum(Q(a;)): 

El = i, ^2 = Sum(l, -), and 

X X 

A = Sum(l, - f Sum(l, - (Sum(l, -f ® Sum(l, 4^))) © Sum(l, Sum(l, -))] ). 

X \ X X X X I 

Finally, ev and c) are extended from X to ev' : ProdSum(X) x N ^ ProdSum(X) 
with [x, n) i-^ x{n) and 5' : ProdSum(X) N as follows. 

(1) For / G X we set D'(/) := d{f) and ev'(/, k) := ev(/, k). 

(2) For f,ge ProdSum(X) we set d'{f(Sg) = d'{f(^g) := max(0'(/), 0'(5), 

ev'(/ © g, k) := ev'(/, k) + ev'{g, k) and ev'(/ g, k) := ev'(/, k) ev'ig, k); 
here the operations on the right hand side are from the field K. 
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(3) For r G N, / e ProdSum(X) define 0'(Sum(r, /)) ^ f)'(Prod(r, /)) := + 1, 

k k 

ev'(Sum(r,/),fc) = ^ev'(/,z) and^ ev'(Prod(r, /), fc) = \{eY'{f,i). 

i=r i=r 

Since ev' and ev, resp. and O', agree on X, we do not distinguish them anymore. 
Subsequently, (ProdSum(X), ev, O) (resp. (Sum(X), ev, O)) is called the product-sum 
sequence domain over X (resp. sum sequence domain over X). 

Example 2.4. The expressions from Example 2.3 are evaluated as 
ev{Euk) = Ei{k) = i, ev(£;2,A:) = £;2(A:) = ^ev(-,i) 

i=i 1=1 



and 



ev(A,fc)=A(fc)=^i^ 



k I ' I 



1=1 



r=l 

We have D{Ei) = 1, (:){E2) = 2 and d{A) = 4. 

Usually, we stick to the following more convenient and frequently used notation. 

• We write, e.g., E = a® Sum(l, cSum(2, 6)) G Sum(X) with a,6, c S X in the 
form 

n i 

E' = cv(a, n) © ^ cv(c, i) ^ cv(6, j) 

i=l j=2 

for a symbolic variable n. Clearly, fixing the variable n, the two encodings E and 
E' can be transformed into each other; if we want to emphasize the dependence 
on n, we also write E' £ Sum„(X). 

• Even more, by abuse of notation, we use instead of © and the usual field 
operations in K. This "sloppy" notation immediately produces the evaluation 
mechanism: cv{E, k) = E(k) for a concrete integer fc G N is produced by substi- 
tution in E' the variable n with the concrete value fc G N. 

• Finally, whenever possible, the evaluation ev(a, n) for some a £ X is expressed by 
well known functions, like, e.g., ev(l/ {x + 1), n) = (n+i) i^'^- 2-1) or ev{xi, n) = 
qf for 1 < i < e (Ex. 2.2). 

Example 2.5. We write Ei,E2,A G I]„(Q(a-)) from Example 2.4 in the more 
convenient notation Ei ^ ^, E2 ~ Hn and (1.1); note that for Ei we must require 
that n is only evaluated for n > 1. 

Let {X, ev, O) be a measured sequence domain and consider the sum sequence 
domain (Sum(X), ev, O) over X. We define the Sum{X)- optimal depth of A G 
Sum(X) as 

min{c)(B)|B G Sum(X) s.t. (1.3) for some A G N}. 
Then we are interested in the following problem. 



DOS: Depth Optimal Simplification. Given A G Sum(X); find B G Sum(X) and 
A G N s.t. (1.3) and s.t. ©(S) is the Sum(X)-optimal depth of A. 



^Note that ev'(Prod(r, /), fc) might be if r is too small. Later, products will be used only as 
described in Ex. 4.10 or in the General case 4.11; there we will take care of the bound r by (4.9). 
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Example 2.6. Consider, e.g., A e Sum(Q(.T)) from Example 2.5. Then with 
our algorithms, see Example 3.13, we find B e Sum(K(.T)) with (1.2) such that 
A{n) = B{n) for all n G N. At this point it is easy to sec that B cannot be 
expressed with dcpth< 1, and thus B \s & solution of DOS. In Section 5 we will 
show that this fact is an immediate consequence of our algebraic construction. 
Summarizing, 2 is the Sum(Q(.T))-optimal depth of A. 

We shall solve problem DOS algorithmically, if X is, e.g., the rational sequence 
domain (Ex. 2.1) or the g-mixed sequence domain (Ex. 2.2). More generally, X 
might be a sequence domain in which a finite number of objects from ProdSum(X') 
can be represented; in this setting X' might be the rational, g-rational or g-mixed 
sequence domain. Note that the General case 4.11 (page 13) includes most of the 
((?-)hypergeometric or g- mixed hypergeometric terms (see Ex. 4.10). 

3. Step I: Reducing the problem to difference fields by telescoping 

Let (X, ev, O) be a measured sequence domain. Then for / G ProdSum(X) and 
r G N the sum S = X]fc=r fi^) ^ ProdSum„(X) satisfies the recurrence relation 

(3.1) S{n+l) = S{n) + f{n + l) Vn > r, 

and the product P = J2'k=r /(^) ^ ProdSum„(X) satisfies the recurrence relation 

(3.2) P{n + 1) = f{n + l)P{n) \/n>r. 

As a consequence, we can define a shift operator acting on the expressions S{n) and 
P{n). Subsequently, we shall restrict to sequence domains X such that the sums 
and products S{n) and P{N) can be modeled in difference rings. 

In general, a difference ring (resp. difference field) (A, a) is defined as a ring A 
(resp. field) with a ring automorphism (resp. field automorphism) cr : A — > A. The 
set of constants consto-A = {fc G A | (j{k) = k} forms a subring'^ (resp. subfield) of 
A. We call constcrA the constant field of (A, a). A difference ring (resp. difference 
field) (E, a) is a difference ring extension (resp. difference field extension) of a 
difference ring (resp. difference field) (A, cr') if A is a subring (resp. subfield) of E 
and (j'(/) = cr(/) for all / G A; we call (A, a') also a sub-difference ring (resp. field) 
of (E, cr). Since a and a' agree on A, we do not distinguish them anymore. 

Example 3.1. For the rational function field K(a;) we can define uniquely the 
automorphism a : K(x) —^ K(x) s.t. cr(x) = x + 1 and s.t. cr(c) = c for all c G K; 
(IK(2;), cr) is called the rational difference field over K. 

Example 3.2. For the rational function field E := K{x, xi,. . . , x,„) from Ex. 2.2 
we can define uniquely the field automorphism cr : F ^ F such that cr(x) ~ x + 1 
and cr(xi) = qtXi for all 1 < i < m and such that cr(c) = c for all c G IK. The 
difference field (E, a) is also called the q-mixed difference field over K. 

Then any expression in ProdSum(K(x)) (resp. in ProdSum( K(x, xi, . . . , Xm))) with 
its shift behavior can be modeled by defining a tower of difference field exten- 
sions over (IK(a;),CT) (resp. of (K(a;, . . . , Xm), cr)). Subsequently, we restrict to 
those extensions in which the constants remain unchanged. We confine to IIS*- 
extensions (SchOl) being slightly less general but covering all sums and products 
treated explicitly in Karr's IlS-extensions (Kar85). 



Subsequently, we assume that consto-A is always a field, which we usually denote by K. Note 
that this implies that Q is a subfield of K. 
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Definition 3.3. A difference field extension (F(i), cr) of (F, a) is called a ns*- 
extension if both difference fields share the same field of constants, t is tran- 
scendental over F, and a{t) = t + a for some a G F* or (j{t) = at for some 
a G F*. If cr{t)/t G F (rcsp. cr(t) — t G F), we call the extension also a H-extension 
(resp. Ti* -extension). In short, we say that {¥(ti) . . . (tg), a) is a nE*-extension 
(rcsp. Il-extension, I]*-extcnsion) of (F, a) if the extension is given by a tower 
of nE* -extensions (rcsp. Il-cxtcnsions, E*-extcnsions). Wc call a nE*-extcnsion 
(F(ii) . . . (te),CT) of (F, cr) with cr(i,) = cti t; + f3i generalized d'Alembertian, or in 
short polynomial, if G F* and /3j G ¥[ti, . . . ,tj_i] for all 1 < z < e. A UY.* -field 
(K{ti) . . . (ig), cr) over K is a nE*-cxtension of (K, a) with constant field K. 

Remark 3.4. If (F(ti) . . . (te), cr) is a polynomial nE*-extension of (F, a) then 
it follows that {¥[ti] . . . [te], cr) is a difference ring extension of (F, cr). 

Karr's approach. The following result from (KarSl) tells us how one can 
design a nE*-field for a given product-sum expression. 

Theorem 3.5. Let (F(t),cr) be a difference field extension o/ (F, cr) with cr(t) = 
at + f where a G F* and / G F. Then the following holds. 

(1) (F(t),cr) is a Ti* -extension of (F, cr) iff a = I and there is no g e¥ s.t. 

(3.3) a{g) = g + f. 

(2) (F(t), cr) is a Il-extension of (F, (7)ifft^0,f = and there are no g g¥* and 
m > such that u{g) = a"^g. 

E.g., with Theorem 3.5 it is easy to see that the difference fields from Exam- 
ples 3.1 and 3.2 are nE*-fields over K. 

From the algorithmic point of view we emphasize the following: For a given 
nE*-field (F,ct) and / G F, Karr's summation algorithm (KarSl) can compute a 
solution 5 G F for the telescoping equation (3.3), or it outputs that such a solution 
in F does not exist; for a simplified version see (Sch05c). In this case, we can adjoin 
a new E*-extension which produces by construction a solution for (3.3). 
Summarizing, Karr's algorithm in combination with Theorem 3.5 enables one to 
construct algorithmically a nE*-ficld that encodes the shift behavior of a given 
indefinite nested sum expression. 

Example 3.6. We start with the nE*-ficld (Q(x), cr) over Q with cr(a-) = x+1. 
Now wc consider the sum expressions of A in (1.1), say in the order 

(3.4) 

i—1 i—1 i—1 i—1 

and represent them in terms of E*-extcnsions following Theorem 3.5.1. 

(1) Using, e.g., Gosper's algorithm (Gos78), Karr's algorithm (KarSl) or a sim- 
plified version of it presented in (Sch05c), we check that there is no g G Q{x) 
with (T{g) = g -\- Hence, by Theorem 3.5.1 we adjoin i7„ in form of the 
E*-extension {Q{x){h),a) of (Q(a;), cr) with a{h) = h + note that the shift 
behavior H„+i = Hn + is refiected by the automorphism cr. 

(2) With the algorithms from (KarSl) or (Sch05c) we show that there is no g G 
Q(x)(/i) with (7{g) =g + ^. Thus we take the E*-extension (Q(x)(/i)(s), cr) of 
{Q{x){h),a) with cr(s) — s -\- and express S by s. 
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(3) With the algorithms from above, we find g = 2s — /i^ G Q{x){h){s) with a{g) — 
g + (^-pTjT. and represent* u'i^^ by g. 

(4) There is no 5 e Q(a;)(/i)(s) with (7{g) = g + '2^^', thus we rephrase T as i in 
the S*-extension {q{x){h){s){t) , a) of {Q{x){h){s),a) with a{t) = i + 

(5) There is no 5 G Q{x){h){s){t) s.t. a{g) = g + ^^^p-] thus we represent A with a 
in the E*-ext. {Q{x){h){s){t){a),a) of {Q{x){h){s){t),a) with ^(a) = a + 

Reformulating a as a sum expression (for more details see Section 4) yields 

(3.5) W ^ 

with ^(n) = W{n) for aU 71 £ N. 

We remark that the sums occurring in W pop up only in the numerator. Here 
the following result plays an important role. 

Theorem 3.7 ((Sch09),Thm. 2.7). Let (F(ti) . . . (t^),a) he a polynomial mi* - 
extension of (F, ct); let A = F[ti, . . . ,te]. Then for all g G A, a{g) — g G A iff 
geA. 

Namely, if, e.g., A consists only of sums that occur in the numerator, then by 
solving iteratively the telescoping problem, it is guaranteed that also the telescoping 
solutions will have only sums that occur in the numerators. 

Remark 3.8. Similar to the sum case, there exit algorithms (KarSl) which can 
handle the product case; for details and technical problems we refer to (SchOSb). 
Note that H-cxtensions will occur later only in the frame of General case 4.11. At 
this point one has explicit control how the sequence domain (X, ev,D) for Sum(X) 
is defined. 

A depth-refined approach. The depth of W in (3.5) is reficctcd by the 
nested depth of the underlying difference field constructed in Example 3.6. 

Definition 3.9. Let (F, a) be a nE*-field over K with F := K{ti) . . . (t^) where 
(7{ti) = aiti or a{ti) = ti + Oi for 1 < i < e. The depth function for elements o/F, 
(5k : IF ^ is defined as follows. 

(1) For any g £ K, 6K{g) ■= 0. 

(2) If (5k is defined for (K(ti) . . . (ii_i), cr) with i > 1, we define 5'^{ti) := Sw^o-i) + 1; 
for 5 — 1^ e IK(ti) . . . {ti), with (?i, 32 G ]K[ti, ■ ■ ■ ,ti\ coprime, we define 

^k(3) niax({(5K(ij)|l < j <i and tj occurs in gi or 52} U {0}). 

The extension depth of a nS*-extension (F(xi) . . . {xr),cr) of (F, (t) is defined by 
max((5K(a:i), . . . , 5^{xr),Q)- 

Example 3.10. In Example 3.6 we have 5q{x) ~ 1, (5q(/i) = 2, (5q(s) = 3, 
(5Q(i) = 4, and ^q{a) = 5. 



Note that there is no way to adjoin a S*-extension of the desired type o-[h2) = h2 + 
l/(x + 1)^, since otherwise a(g — /i2) = {g — ^2), i.e., consto-Q(a;)(/i)(s)(/i2) ^ Q. 
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With the approach sketched in Example 3.6 we obtain an alternative sum rep- 
resentation W{n) for A{n) with larger depth. Motivated by such problematic situa- 
tions. Karr's ni]*-fields have been refined in the following way; see (SchOSa; Sch08). 

Definition 3.11. Let (F, a) be a nS]*-field over K. A difference field extension 
(F(s), a) of (F, a) with a{s) = s+f is called depth-optimal I]*-extension, in short S*- 
extension, if there is no S]*-extension (E, a) of (F, a) with extension depth < (5k(/) 
such that there is a g S E as in (3.3). A nS*-extension (F(f;i) . . . (te), cr) of (F, cr) 
is depth-optimal, in short a HS'^-extension, if all S*-extensions are depth-optimal. 
A nS^-field consists of 11- and S'^-extensions. 

Note that a E'^-extension is a I]*-extension by Theorem 3.5.1. Moreover, a nE*-field 
(F,(t) with depth < 2 and x G F such that cr{x) = x -I- 1 is always depth-optimal; 
see (Sch08, Prop. 19). In particular, the rational and the q-mixed difference fields 
from the Examples 3.1 and 3.2 are WYi^ -fields over K. 

Given any HS'^-field, we obtain the following crucial property which will be essential 
to solve problem DOS. 

Theorem 3.12 ((Sch08), Result 3). Let (F,cr) be aUJ:^ -field overK. Then for 
any f,g£¥ such that (3.3) we have 

(3.6) (5k(/) < SKig) < <5k(/) + 1. 

In other words, in a given HS'^-field we can guarantee that the depth of a telescoping 
solution is not bigger than the depth of the sum itself. 

Example 3.13. We consider again the sum expressions in (3.4), but this time 
we use the refined algorithm presented in (SchOS). 

(1) As in Example 3.13 we compute the HS'^-field {Q{x){h), a) and represent _ff„ 
with h. From this point on, our new algorithm works differently. 

(2) Given {Q{x){h),a), we find the E-'-extension {Q{x){h){h2),(T) of {Q{x){h),a) 
with cr(/i2) = /i2 + (^x+iy which we find s' = ^{h^ -f /12) such that cr(s') — s' = 

Hence we represent S by s'. 

(2) X 

(3) Hn can be represented by /12 in the already constructed HE -field. 

(4) Our algorithm finds the E'^-extension (Q(.T)(/i)(/i2)(/i3), cr) of {Q{x){h){h2), cr) 
with cr(/i3) = /i3 + (^,j.^iy.i together with t' = ^{h^ + 3hh2 + 2/13) such that 

cr{t') — t' = hence we rephrase T as t'. 

(5) Finally, we find the E-'-ext. (Q(x)(/i)(/i2)(^3)(/j4), cr) of (Q(a;)(/i)(/i2)(/i3), ct) 
with cr(/i4) = /i4 -I- (j^pYF and get a' = ^{h^ -\- 2h^ -I- 6{h -|- 1)^.2/1 -l- Sh^ -\- {8h -|- 

4)ft,3 + 6/14) s.t. cr(a') — a' = ^^|qrr^; ^ is represented by a'. 
Reinterpreting a' as a sum expression gives B in (1.2); see also Example 2.6. 

To sum up, we can compute step by step a HE'^-field in which we can represent 
nested sum expressions. To be more precise, we will exploit the following 

Theorem 3.14 ((Sch08),Result 1). Let (F, a) be a UE^ -field over K and f e F. 

(1) There is a T,^ -extension (E, cr) of (F, cr) in which we have 5 S E such that (3.3); 
(E, cr) and g can be given explicitly ifK. has the form as stated in Remark 3.15. 

(2) Suppose that (F, cr) with F = (G-(yi, . . . , y,.) is a polynomial WTi^ -extension of 
(G, a) . If f £ G[?/i . . . ,yr\, then (E, cr) from part (1) can be given as a polynomial 
HT,^ -extension of (G, cr); if¥. = F(ii, . . . , te), then g e G[2/i, . . . , yr][ti, . . . ,te]. 
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Remark 3.15. From the computational point of view certain operations must 
be carried out in the constant field'''. For instance, Theorem 3.14 is completely 
constructive, if IK is of the following from: K = A((7i, . . . , g^) is a rational func- 
tion field with variables gi, . . . , over an algebraic number field A. Due to the 
restrictions of the computer algebra system Mathematica, the implementation in 
Sigma (Sch07) works only optimal if A = Q, i.e., K = Q{qi, ■ ■ ■ , Qm)- 

4. Step II: Reinterpretation as product-sum expressions 

Let (X, ev,5), e.g., be the g-mixed sequence domain from Example 2.2 with 
X — K{x, xi, . . . , Xm) and let (F, a) with X = F be the g-mixcd difference field 
from Example 3.2. Moreover, take A G Sum(X). 

Then in the previous section we have demonstrated how one can compute a 
polynomial S'^-extension (F(ti) . . . (tg), (t) of (F, a) in which one can model the shift- 
behavior of A by an clement a G F[ti, . . . , tg]. Then, as illustrated in Example 3.13, 
we were able to reinterpret a as an element from B £ Sum(X) such that (5k (a) = 
and such that (1.3) where A £ N could be given explicitly. 

In order to accomplish this task algorithmically, we will supplement the con- 
struction of the difference field (F(ii) . . .{te),a) by defining in addition an explic- 
itly given difference ring monomorphism. Namely, following (Sch09) we will em- 
bed the difference ring {¥[ti, . . . ,te],cr) into the ring of sequences by a so-called 
K-monomorphism. It turns out that any element h € F[ti, . . . ,te] can be mapped 
injectively to (cv(iJ, k))k>o for some properly chosen expression H € ProdSum(X). 

Subsequently, we define the ring of sequences and K-monomorphisms. Let K be 
a field and consider the set of sequences with elements (a„)„>o = (oo, ai, ^2, • ■ • ), 
fli G K. With component-wise addition and multiplication we obtain a commuta- 
tive ring; the field K can be naturally embedded by identifying fc G K with the 
sequence {k, k, k, . . .); we write = (0, 0, 0, ... ). 

We follow the construction from (PWZ96, Sec. 8.2) in order to turn the shift 
(4.1) S : (ao,ai,a2,...) (ai, a2, 03, . . . ) 

to an automorphism: We define an equivalence relation ^ on with (a„)„>o ~ 
{bn)n>o if there exists a c? > such that Ofc = bk for all k > d. The equivalence 
classes form a ring which is denoted by S'(K); the elements of S{K) (also called 
germs) will be denoted, as above, by sequence notation. Now it is immediate 
that S : S(K) — * S{K) with (4.1) forms a ring automorphism. The difference ring 
{S(K),S) is called the ring of sequences (overM.). 

A difference ring homomorphism r : Ai ^ A2 between difference rings (Ai, ai) 
and (A2, 0-2) is a ring homomorphism such that r(cri(/)) = (T2(r(/)) for all / G Ai. 
If T is injcctive, we call r a difference ring monomorphism. 

Let (A, cr) be a difference ring with constant field K. Then a difference ring 
homomorphism (resp. difference ring monomorphism) r : A — > iS'(K) is called K- 
homomorphism (resp. ¥^-monomorphism or 'K-embedding) if for all c G K we have 
that r(c) = (c, c, . . . ). 

As mentioned already above, our final goal is to construct a K- monomorphism 
T : F[ti, . . . , te] S{K). For this task we exploit the following property. 



Actually, we require the same computational properties as for Karr's summation algorithm; 
see (KarSl, Thm. 9) or (KS06, Thm. 3). 
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If r : A — > S{K) is a K-homomorphism, there is a map ev : A x N — > K with 

(4.2) r(/) = (ev(/,0),cv(/,l),...) 

for all / e A which has the following properties: For all c G K there is a 5 > with 

(4.3) yi>6:ev{c,i) = c; 

for all /, g G A there is a 5 > with 

(4.4) Vi > (5 : ev(/ g, i) = ev(/, i) evig, i), 

(4.5) Vi > 5 : ev(/ + g,i) ^ ev(/, i) + ev{g, i); 

and for all / G A and j G Z there is a (5 > with 

(4.6) yi>6: ev(cr^(/),0 = cv{f,i + j). 

Conversely, if there is a function ev : A x N ^ IK with (4.3), (4.4), (4.5) and (4.6), 
then the function r : A — > S(K) defined by (4.2) forms a K-homomorphism. 

Subsequently, we assume that a K-homomorphism/K-monomorphism is always 
defined by such a function ev; ev is also called a defining function of r. To take 
into account the constructive aspects, we introduce the following functions for ev. 

Definition 4.1. Let (A,(t) be a difference ring and let r : A ^ S{K) be a K- 
homomorphism given by the defining function ev as in (4.2). ev is called operation- 
bounded by L : A ^ N if for aU / G A and j G Z with S = S{f,j) L{f) + 
max(0, — j) we have (4.6) and for all /, g G A with 5 = S{f,g) :~ niax{L{f), L{g)) 
we have (4.4) and (4.5); moreover, we require that for all / G A and all j G Z we 
have L{a^{f)) < L{f) + max(0, —j). Such a function is also called o-function for 
ev. ev is called zero-hounded by Z : A — > N if for all / G A* and all i > Z{f) we 
have ev(/, i) ^ 0; such a function is also called z-function for ev. 

Example 4.2. Given the HE-'-field (K(x),cr) over K with (7{x) = a; + 1, we 
obtain a K-homomorphism r : IK(a;) — > S{K) by taking the defining function (2.1); 
here we assume that / = ^ G ]K(a;) with p G K[a:;] and q G K[a:]* are coprime. For 
the o-function L{f) we take the minimal non-negative integer I such that q{k-\-l) ^ 
for all A: G N, and as ^-function we take ^(/) = L{pq). Note: Since p{x) and q{x) 
have only finitely many roots, t(^) = iff ^ = 0. Hence r is injective. 

bummarizmg, the nS^-field (IK(a;), cr) with a{x) = x 1 can be embedded into 
{S{K),S). More generally, if (F, ct) is the q-mixed difference field, t : F -> S{K) 
with the defining function ev given in (2.2) is a K-monomorphism. In addition, 
there are a computable o-function L and a computable z-function Z for ev; for 
details we refer to (Sch09, Cor. 4.10) which rehes on (BP99). 

Example 4.3. Take the rational difference field (K{x),a) and the K-monomor- 
phism r with defining function ev and the o-function L from Example 4.2 and con- 
sider the E"'-extension {K{x){h),a) of {K{x),a) with a{h) = h + We get a 
K-homomorphism r' : K(x)[/i] —>■ S{K) where the defining function ev' is given by 
ev'(/i, k) = Hk and 

d d 

ev'(^ f^h\ k) = J2 ev(/., k)ey'{h, k)\ 

i=0 i=0 
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As 0-function we can take L' {^'^^^^ fih"^) = max(L(/i)|0 < i < d). Now suppose 
that r' is not injective. Then wc can take / = SiLo/«^* ^ IK(a;)[/i] \ {0} with 
deg(/) = d minimal such that r'(/) = 0. Since r is injective, / ^ K(a;). Define 
g := <j{fd)f — fd<^{f) G K{x)[h]. Note that deg{g) < d by construction. Moreover, 

r'(5) = rK/,))r'(/)-r(/,)r'K/)). 

Since t'(/) = by assumption and r'(cr(/)) = 5(r'(/)) = 5(0) = 0, it foUows 
T'{g) = 0. By the minimaHty of deg(/), g = 0, i.e., a{fd)f - fd(j{f) = 0, or equiv- 
alcntly, ^ G ]K(a;). As / ^ IK(a;), this contradicts (KarSl, Theorem 4). 

Example 4.4. Take the rational difference field (K(a;), cr) and the K-monomor- 
phism T with defining function ev and o-function L from Example 4.2, and con- 
sider the n-extcnsion (K{x){b),a) of (K{x),a) with cr(fe) = 2(2^+1) ^- 
a K-homomorphism t' : ]K(a;)[6] ^ S'(K) with its defining function specified by 
ev'(b,k) = 11^=12(2^1) ^ CkV^ ^^'^^ (4-^) where t := b; note that r'(6) has 
no zero entries by construction. We take -^'(X^iLo /i^*) = inax(i(/i)|0 < i < d)) as 
o-function. By similar arguments as in Ex. 4.3 it follows that r is injective. 

More generally, we arrive at the following result; see (Sch09) for a detailed proof. 

Lemma 4.5. Let (F(ii) . . . {te){t),(7) be a polynomial HY,* -extension of (F, ct) 
with K := const^F and a{t) ^ at + p. Let t : ¥[ti] . . . [te] S{K) he a K- 
monomorphism with a defining function ev as in (4.2); let L be an o-function for 
ev and let Z be z-function for ev\f (ev is restricted on ¥). Then: 

(1) There is a M.-monomorphism t' : F[ti] . . . [te][t] — > S{K) with a defining function 
ev' such that ev'|F[(j = ev; if P ~ 0, ev'(t,k) ^ for all k > r for some 
r S N. Such a r' is uniquely determined by 



(4.7) ev'{t,k) 



k 

C 



J^ev(a,z — 1) if a{t) — at 



i—r 
k 



'^ev{P,i-l) + c ifa{t)=t + P, 



i—r 



up to the choice 0/ r G N and c G IK; we require c 7^ 0, if P ~ 0. 

(2) Fixing (4.7) we obtain, e.g., the following defining function for t' : 

d d 

(4.8) ev'(^/.f,fc) :=^ev(/„fc)ev'(t,fc)' VA; G N. 

(3) In particular, there is an o-function L' for ev' with -Zj'|F[ti,...,te] — '"^'^ ^ 
are computable, L' can he computed. We can choose (as a constructive example ) 



(4.9) 



max(L(Q!), Z(q;)) + 1 if a{t) ^ at 
L(/3)-fl ifa{t)^t + p. 



Remark 4.6. Let (F(ti) . . . (te), f) be a polynomial ns*-ext. of (F, cr) with 
K = constcrF and let r : F[ti] . . . [te] — > S{K) be a K-homomorphism with a defining 
function ev. Then there is implicitly a z-function for evjip-; see (Sch09, Lemma 4.3). 

Applying Lemma 4.5 iteratively produces the following result. 
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Theorem 4.7. Let (F(yi) . . . {yr){ti) . . . (te), cr) be a polynomial HE* -extension 
of (F,cr) with K := const^F; let t:A^S{K) with A = ¥[yi,...,yr] be a K- 
embedding with a defining function (4.2) and with an o-function L. 
Then there is a "K-embedding r' : h.[ti\ . . . [te] — > S'(K) with a defining function ev' 
and with an o-function L' s.t. cv'|a ~ cv and L'|a = L. 

This construction that leads in Theorem 4.7 to ev' is called canonical, if it is per- 
formed iteratively as described in (4.7) and (4.8) of Lemma 4.5. Note that any 
defining function ev' with (4.7) evaluates as (4.8) if k is chosen big enough; in our 
canonical construction we assume that (4.8) holds for all fc > 0. 

Since r' in Lemma 4.5 is uniquely determined by (4.7) the following holds. 

Theorem 4.8. Let (F, a) be a WE* -field over K and let (F(?/i) . . . {yr),o) be a 
polynomial HE* -extension of (F, a). Let r : F[yi, . . . , y^] ^ >5'(K) be a IK-embedding 
with a defining function ev together with an o-function for ev. Take the measured 
sequence domain (F, ev, i5k). Then for any f G ¥[yi, . . . ,yr], there is an F £ 
ProdSum{¥) s.t. T{f) = {F{k))k>o and 0(F) = 5k(-F). 

Summarizing, given such a difference ring (F[yi, ...,?/,,], cr) and IK-monomorphism, 
one can rephrase the elements of F[yi, . . . ,yr] as expressions from ProdSum(X) 
such that the depth of both domains are identical. 

Example 4.9. Take the nS'^-field (Q(x),(t) with a{x) = .t + 1 together with 
the Q-embedding r : Q{x) S{Q) with defining function (2.1) as carried out in 
Ex. 4.2 (K = Q); let (Sum(Q(a;)), ev, O) be the sum sequence domain over Q(x). 
Moreover, consider the nS-^-field (Q(a;)(/i)(/i2)(/i3)(/i4), cr) from Ex. 3.13. Then we 
can construct the Q-embedding r' : Q(x)[/i, /12, ^ S{Q) with the defining 
function ev' which is canonically given by cv'Iq^^.) = ev and by 

ev'(/i,fc) =cv(Sum(l,i),fc) =i7fc, 

(4.10) 

eY'{hj,k) ^ ev{Sm-n{l,^),k) ^ Hlf^ for j 6 {2,3,4}. 

Note that t)(Sum(l, i)) = (5q(/i) and 5(Sum(l, i)) = 5q(/ij) for j G {2,3,4}. 

Recall that we want to solve problem DOS for a measured sequence domain 
{X, ev, D). In Section 5 we shall solve this problem for the following general setting. 
Loosely speaking, the terms of indefinite nested sums and products X are modeled 
by polynomials from F[?/i, . . . , i/r] and the reinterprctation of the corresponding 
product-sum expressions is accomplished by its K-monomorphism from F[?/i, . . . , y^] 
into the ring of sequences; in particular, the depth of such a product-sum expression 
is equal to the depth of the corresponding polynomial from F[?;i, . . . , y^]- 

Example 4.10. We start as in Example 4.9, but now we take the H-extension 
(F(yi) . . . {yr),<7) of (F, cr) such that = e F for 1 < j < r. Then we can 

extend the K-embedding r to r' : F[yi, . . . ,yr] S{K) with the defining function 
ev' canonically given by ev'jp = ev and 

fc 

ey'{yj,k) = Cj J| a.j{i) 

i—Tj 

for all 1 < .7 < 7" with rj > Z{aj) and Cj G K*. Note: with Fj = CjProd(rj, aj) we 
have Fj{k) = ev'{yj,k) and t)(Fj) = SK^yj)- Moreover, we can model a finite set of 
hypergeometric terms in the sequence domain {X, ev', Sk) with X :=¥[yi, . . . ,yr]. 
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Similarly, wc are in the position to handle g-hypergeometric sequences or mixed 
hypergeometric sequences. More generally, we can handle the following case. 

GENERAL CASE 4.11. The ground field". Let (F, a) be a nS-^-field over K, 
let To : F — *■ S(K) be a K-embcdding with a defining function evo : F x N ^ K, 
and let Lq : 1^ — * N be an o- function and Z : F ^ N be a z-function for evo; 
moreover, consider the sequence domain (F,evo,i5K). 

A polynomial extension. In addition, choose a polynomial HS'^-cxtcnsion 
(F(yi) . . . (?/,.),cr) of (F,cr) and set X := F[?/i, . . . , y^]. Then extend the K- 
embedding tq to t : X ^ S(K) by extending the defining function cvq canon- 
ically to ev : X x N ^ K and by extending the o-function Lq to L following 
Lemma 4.5; if Lq is computable, also L is computable. By construction it fol- 
lows that for 1 < i < 7' there exist Fi G ProdSum(F) such that 

(4.11) ev{y,,k)^F,{k) Vfc > and ^F,) = SK{y^)■ 

In particular, for each f G X, one gets explicitly an F G ProdSum(X) s.t. 
ev(/, k) = F{k) for aU fc G N and s.t. D{F) = Smif). 

The sequence domain. We obtain the sequence domain {X, ev, (5k) which 
models the product-sum expressions (4.11) with the depth given by S^. 

"E.g., wo can take the g-mixed difference field (F, cr) with F = ¥^(x, xi, . . . , Xm) from Ex. 3.2, 
and we can take the K-embedding tq : F — > S(K) where evg := ev is defined as in (2.2); note 
that the measured sequence domain (F, evo, (5k) been presented in Ex. 2.2. From the 
computational point of view we assume that K is of the form as stated in Remark 3.15 



5. Combining the steps: Finding optimal nested sum representations 

E.g., for the g-mixed sequence domain (F, ev,()) from Ex. 2.2 with X = F = 
K{x,xi, . . . ,Xm) we will solve problem DOS for A G Sum(F) as follows; here we 
assume that K is of the form as stated in Remark 3.15 

Take the 5-mixed difference field (F, cr) over K with the automorphism a de- 
fined in Example 3.2. Moreover, take the the K-embedding r : F — > S{K) with the 
defining function ev given in (2.2), and choose a computable o-function L and a 
computable z-function Z for ev. Then by Theorem 5.1 below the following con- 
struction can be carried out algorithmically. 

Step I: Reduction to a LtS^-field. Given A G Sum(F), construct a polyno- 
mial E'^-extcnsion (F(.si) . . . (s„),(t) of (F, cr) and extend the K-monomorphism r 
to t' : F[si, . . . , Sti] S{K) with a defining function ev' such that the following 
additional property holds: We can take explicitly an a G F[si, . . . , s„] and a A G N 
such that 

(5.1) ev'(a, k) = A{k) Vfc > A and^ SK{a) < d{A). 

Note that we rely on the fact that all our sums are represented in HS'^-fields; for 
general nS*-fields 6-K{a) might be bigger than <}{A), see Example 3.6. 

Step II: Reinterpretation as a product-sum expression. In particular, by the 
concrete construction of the K-monomorphism based on the iterative application 
of Lemma 4.5, construct a. B E Sum(F) such that 

(5.2) cv'{a,k) ^ B{k) Vfc > and 5k (a) = 5 (S). 



^Recall that A(k) = ev(j4, fc); here ev is the evaluation function of the sequence domain 
{X, ev, 0) where X = ¥ (or X = V[yi, . . . , yr] as defined in the General case 4.11). 
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Then due to the properties of the HS'^-field and the fact that r' is a IK-mono- 
niorphisni (in particular, that t' is injective), we wiU show in Theorem 5.5 that the 
depth of _B G Sum(F) is Sum(F)-optimal. i.e., B together with A are a solution of 
problem DOS. 

We will solve problem DOS for the General case 4.11 by applying exactly the 
same mechanism as sketched above. 

Theorem 5.1. Let (F(yi) . . . (y,.), ct) with X = ¥[yi,...,yr] be a UT.^ -field 
over K, let t : X ^ S'(K) be 'K- embedding with ev, L and Z , and let {X, ev, (5k) be a 
sequence domain as stated in the General case ^-11; in particular let (Sum{X), ev, d) 
be the sum sequence domain over X . Then for any A £ Sum{X) there is a S*- 
extension D := (F(yi) . . . (yr){si) . . . (su), <j) of (¥(yi) . . . (j/r), cr), where D is a 
polynomial extension of {¥,a), and there is a K.-embedding 

(5.3) t':X[.si,...,s„]^5(K) 

where the defining function ev' and its o-function L' are extended from X to 
X[si, . . . , Sy], with the following property: There are A G N and a £ A such 
that (5.1); in particular, for any h G ^[si, . . . , s^] there is an H G Sum{X) s.t. 

(5.4) ev'(/i, k) ^ H{k) Vfc > and Sw^h) = d{H). 

This extension, the defining function ev' for r', A, and a can be given explicitly, if 
L and Z are computable and i/K has the form as stated in Remark 3.15. 

Proof. We show the theorem by induction on the depth. If A G Sum(X) 
with X){A) = 0, then A G IK and the statement clearly holds. Now suppose that we 
have shown the statement for expressions with depth< d and take A G Sum(X) 
with X){A) = d+\. Let Ai^ . . . ,Ai be exactly those subexpressions of A which do 
not occur inside of a sum and which cannot be split further by © and ®, i.e., for 
1 < i <l, either G X or Ai is a sum. First we consider Ai. If Ai G X, then for 
ri = 0, ai = j4i, ev' = ev and i = 1 we have 

(5.5) ev'(aj, k) = A^{k) Vfc > r^, and 5K{at) < ^{A^). 

Moreover, the property (5.4) for any h £ X holds by choosing H := h. 
Otherwise, Ai = Sum(Ai,i^i) for some Ai G N and Fi G Sum(X) with d{Fi) < d. 
If 0(^i) < d, we get by induction a E'^-extcnsion D :— {¥{yi) . . . {yr){si) . . . (su), cr) 
of (F(?/i) . . . {yr),(y) where D is a polynomial extension of (F,cr); moreover we can 
extend r to a K-embedding t' : D — + S'(K) with n := X[si, . . . , Su] where its defin- 
ing function ev' is extended from X to D, and we can extend the o-function L to an 
o-function L' for ev' s.t. the following holds: There are oi G B and ri G N s.t. for 
i = 1 we have (5.5); in particular, for any /i G D, there is an iJ G Sum(X) s.t. (5.4). 
If 5(^i) = d + 1, we can take by the same reasoning such an extension D of 
(F, cr) with D := X[si, . . . , Su] and r' with a defining function ev' together with an 
o-function L' in which we can take /i G B and /i G N s.t. 

ev'(/i,fc) =Fi(fc) Vfc>/i, and,5K(/i) <c)(Fi). 

By Theorem 3.14 take a S'^-extension E := (F(yi) . . . {yr){si) . ■ . {su){ti) . . . (it.),cr) 
of (F(j/i) . . . (yr)(si) . . . (su), cr) such that E a polynomial extension of (F, a) and in 
which we have 5 G E with E := D[ii, . . . ,ty] such that 

(^{9) =5 + o'(/i)- 
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Moreover, by iterative application of Lemma 4.5 we can extend the K-cmbcdding 
t' from D to a K-embedding r' : E ^ S{K) by extending the defining function ev' 
canonically from D to E, and we can extend L' to an o- function for ev'; note that 
this construction can be performed such that for any /i G E there is _ff G Sum(X) 
with (5.4). Now take^ n := max(/i, L'(/i), + 1) and define c := J2l'Jh -^('^) " 
ev'{g, 7'i — 1) G K. Then for aU ??. > ri, 

ev'{g + c, n) = cv{a^^{g) + fi,n) + c = ev'{g,n- 1) + ev'(/i,7i) + c 

n 

= ev'(5, n-l)+Fi{n) + c^ ■■■ ^ ev\g, n - 1) + ^ Fi(/c) + c Ai{n). 

k—ri 

Set oi := .9 + c G E. Then (5.5) for i = 1. Moreover, since SK{g) < i5k(o'(/i)) + 1 
by Thm. 3.12, 

SKiai) = 5K{g) < <5k((t(/i)) + 1 = (5k(/i) + 1 < f (i^i) + 1 < + 1. 

We continue to consider A2, . . . , Ai and finaUy arrive at a polynomial E'^-extension 
(F(yi) . . . (?/r)(si) . . . {su){ti) ■ ■ ■ (<!,'), cr) of (F, (t); during this construction we can 
extend t to a K-embedding t' : A ^ S(K) with A := X[si, . . . , s„][ti, . . . , tyi] and 
with a defining function ev' and we can extend L to an o-function U for ev' s.t. 
the following holds. We can take ai, . . . , a; G A and ri, . . . , n G N s.t. for 1 < i < I 
we have (5.5); in particular, for any /i G A, there is an iJ G Sum(X) s.t. (5.4). 
Finally, we construct a G A by applying in the given A the substitution Ai ai 
for all 1 < z < Z and by replacing ® and © with the field operations • and +, 
respectively. Then it follows that (5.1) for A := max(ri, . . . , r/, L'(ai), . . . , L'{ai)). 
This completes the induction step. Note that all the construction steps can be 
carried out by algorithms if L and Z are computable and if K has the form as 
stated in Remark 3.15. In particular, ev', a and A can be given explicitly. □ 

Example 5.2. For the input sum (1.1) the presented procedure in Theorem 5.1 
carries out simultaneously the constructions from Examples 3.13 and 4.9: We ob- 
tain the nS-^-field {Q{x){h){h2){hz){hA) , a) over Q together with the Q-embedding 
r' : Q(x)[/i, /i2, /i3, /i4] — *■ 5'(Q) where the defining function ev' is canonically given 
by ev'|Q(j;) = ev for ev as in (2.1) and by (4.10). By construction, we can link the 
sums in (3.4) with h, s', /i2, t', a' G Q(a;)[/i, ^3, ^4] from Example 3.13: for n G N, 

ev'{h,n)~Hn, cv' {s' , n) = S {n) , 

^^■^^ cv'{h2,n)^ H^^\ ev'(i',n) =T(n), cv' {a' , n) ^ A{n) . 

(2) 

The depths of Hk,S,Hi, ,T,A are 2,3,2,3,4, respectively. The corresponding 
depths (5q(/i) = ^q(s') = i5q(^2) = ^q(^') = ^Q(a') = 2 in the HS^-field are the 
same or have been improved. Using (4.10) we can reinterpret, e.g, s' as the sum 
expression F := i(Sum(l, l/a:)^ + Sum(l, l/cc^)) g Sum(Q(x)) with F{n) = S{n) 

for aU n > and d{F) = Sq{s'); this leads to the identity s{n) = ^(i/^ + H^'^). In 
the same way we obtain sum expressions in Sum(Q(a;)) for the IlS'^-field elements 



Remark: If /i is a hypergeometric term, ri can be also obtained by analyzing only /i 
(without knowing the telescoping solution g); for more details see (AP05). 
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h, h2, t\ a' in (5.6), and we arrive at the following identities. For n > 0, 

H^?^=H^?\ r(n) = i(H„3 + 3i/(^)i/„ + 2H(3)), Ain)^Bin) 
where B is given as in (1.2). 

We remark that this translation mechanism presented in Theorem 5.1 is imple- 
mented in the summation package Sigma. Namely, given the General case 4.11 (K 
as stated in Remark 3.15) and given® A G Sum(X), Sigma computes the following 
ingredients: 

• A E'^-extension (F(yi) . . . (yr)(si) ■ • ■ (sti),tT) of (F(?/i) . . . {yr),a) and a K-em- 
bedding (5.3) with a defining function ev' and o-function L' . 

• a £ X[si, . . . , s„] and A £ N such that (5.1). 

• T'(a) is given explicitly by a^ _B G Sum(X) s.t. (5.2). 

Then Sigma outputs for a given A G Sum(X) the result B G Sum(X) with A. 

The final goal is to prove Theorem 5.5 which guarantees that the output B is 
indeed a solution of problem DOS. We start with the following 

Lemma 5.3. Let (F(yi) . . . (y^), ct) be a UY.^ -field over K and let t be a K.- 
monomorphism with ev and L as in the General case ^ei (F(fi) . . . (tg), a) be a 

polynomial I]* -extension of (F, a) and let p : F[<i, . . . , tg] —>■ S{K) be a 'K-embedding 
with a defining function evp and an o-function Lp such that evpl^ = ev. 
Then there is a T,* -extension D := (F(?/i) . . . (?/r)(-Zi) • ■ . {zi),cf) of (F(yi) . . . (i/r), ct) 
where D is a polynomial Y,* -extension of (F, a) with the following properties: 

(1) There is a difference field monomor. (j) : ¥{ti, . . . ,te) —>■ F(yi, . . . , yr){zi, . . . , Zi) 
such that for all a G F(ii, . . . , t^), 

(5.7) SMa)) < S^ia) 
and such that for all a G F[fi, . . . ,t^], 

(5.8) (/)(a) e¥[yi,...,yr][zi,...,zi]. 

(2) There is a M^- embed, t' :¥[yi, . . . ,yr][zi, . . . , zi] S'(K) with a defining function 
ev' and an o-function L' s.t. ev'|j-[j,^ j,^] = ev and s.t. for all a G F[ti, . . . , tg,], 

(5.9) T'(0(a))=p(a). 

Proof. The base case e = holds with (p{a) = a for all a G F and r' := p. Sup- 
pose the lemma holds for e extensions (H, cr) with H = ¥{ti) . . .{tg) and let (ED, a) 
with D := ¥{yi) . . . {yr){zi) . . . (z;), r' with ev' and L', p with evp, and (j) as stated 
above; set ¥ = ¥[yi, . . . ,yr, zi, . . . , zi]. Now let {M{t),a) be a S*-ext. of (H, a) with 
/ := a(t) - t G ¥[ti . . . ,te], and take a K-embedding p' : F[ii, . . . ,ie][i] ^ ^(IK) 



In Sigma A is inserted, e.g., in the form (1.1) without using evaluation functions like (2.1) 
or (2.2); this implies that the lower bounds of the involved sums and products must be chosen in 
such a way that no zeros occur in the denominators during any evaluation. 

^In Sigma the lower bounds of the sums and products are computed by (4.9). Looking 
closer at this construction, no zeros occur in the involved denominators of B when performing the 
evaluation B{n) for n > A. Hence the output can be returned, e.g., in the form like (1.1) or (1.2) 
which is free of any explicit evaluation functions like (2.1) or (2.2). 
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with a defining function p' and an o-function Lpi s.t. evp'|F[(j_ = evp. 
Case 1: If there is no g G D such that 

(5.10) = (/)(/), 

we can take the I]*-cxtcnsion (©(y),^) of (D, cr) with a{y) = y + 4>{f) by The- 
orem 3.5.1 and we can define a difference field monomorphism cj)' : M(t) — > 
s.t. 4>'{a) = 4){a) for aU a G H and such that 4i'{t) = y. By construction, 
5K{y) = SKiHf)) + 1- Since 

(5.11) SKWf)) + l<SK{f) + l = 5Kit), 

6-K{(f){a)) < 6^(0.) for aU a G M{t). Moreover, since (j){f) G E, it follows that 
(E{y),a) is a polynomial extension of (F, cr). Moreover, for all a G F[fi, . . . ,te,t], 
(j){a) G E[j/]. This proves part (1). 

Now we extend the K-embedding r' from E to r' : K[y] S{K) with the defining 
function ev' where ev'{f,k) = ev{f,k) for all / G E and ev'{y,k) is defined as in 
the right hand side of (4.7); here we choose (3 = (j){f) and c = evp'{t,r — 1) for 
some r G N properly chosen. In particular, we can extend the o-function L' for our 
extended ev' by Lemma 4.5 (note that there is a z function for ev' restricted on F 
by Remark 4.6). Then for all fc > r {r is chosen big enough with L'p and L'), 

k 

ev'(0'(t), fc) = ev'(y, fc) - J] ev'(^(/), i - 1) + cvp> (t, r ~ 1) 

(5-12) 

= ^evp(/,i-l)+evp,(t,r-l)= ^ evp{f,i-l) + h{r) 

i—r i—r-\-l 

with h{r) = evp(/, r - 1) + evp- (t, r - 1) = cvp/ (/ 4- r - 1) = cvp/ (cr(t), r - 1) = 
evp/(t,r). Applying this reduction fc — r -I- 1 times shows that 

k 

ev'((/)'(t),fc) = J2 cvp(/,i-l)+evp,(i,r) = ...=evp,(i,fc). 

i=r+l 

Hence t'(0'(<)) = p'(i), and thus T'{(f>'{a)) = p'(a) for all a G F[ti . . . ,te][t]. 

Case 2: Otherwise, if there is a 5 G D s.t. (5.10), then g G E by Theorem 3.7. In 

particular, (5k(5) < (5k('/)(/)) + 1 by Theorem 3.12. With (5.11), it follows that 

(5.13) dK{g) < SKit). 

Now observe that 

S{p'{t)) - p'{t) = p'{a{t) -t) = p{f) = r'(0(/)) = r'{<j{g) - g) = 5(r'(g)) - r'{g). 

Hence S {p' (t) - t' (g)) = p'{t)-T'{g), i.e., p'{t)-T'{g) is a constant in S'(K). Thus, 
p'(i) = T'{g) + (c}„>o for some c G K. Since {(l){M){g), a) is a difference field (it 
is a sub-difference field of (D, fi)), g is transcendental over 0(H) by Theorem 3.5.1. 
In particular, we can define the difference field monomorphism 0' : M{t) — > D with 
(j)'{a) = (j){a) for all a G H and 4)'{t) = g + c. Since g G E, 4)'{t) G E, and therefore 
(j){a) G E for all a G F[ii, . . . , ig] [t]. With (5.13) and our induction assumption 
it follows that SK{T'{a)) < SK{a) for all a G M{t). This proofs part (1). Note 
that T'{(l)'{t)) = T'{g) + c = p'{t) by construction. Hence, T'{(j)'{a)) = p{a) for all 
a G F[ti . . . ,te][i]. This proves part (2) and completes the induction step. □ 
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Note that the proof of Lemma 5.3 is constructive, if the miderlying o- and z- 
functions of r are computable and if K is given as in Remark 3.15. For simphcity, 
ingredients hke r\ (needed, e.g., for (5.12)) have not been specified exphcitly. 

In Theorem 5.4 we can show the fohowing: The nested depth of an element in 
a nS'^-field over K is smaller or equal to the depth of an element in a nE*-field 
over K provided that both elements can be mapped to the same sequence s G S'(K) 
by appropriate K-monomorphisms. 

Theorem 5.4. Let (F(?/i) . . . (y^), cr) he a UT,^ -field over IK and let t be a K- 
embedding with ev and L as in the General case Let (F(ii) . . . (ie),(T) be a 

polynomial Y,* -extension of (F, cr) and let p : F[ti, . . . , tg] — > S'(K) be a 'K-embed- 
ding with a defining function evp and an o-function for evp. Then for any s G 
t(F[2/i, ...,yr])n p(F[ti, . . . , te]) we have 

Proof. Take a I]*-cxtension {¥{yi) . . . {yr){zi) . . . {zi),a) of (F(t/i) . . . {yr),cr), 
a monomorphism : F(ii) . . . {te) ¥{yi) ■ ■ ■ {yr){zi) ■ ■ ■ {zi) s.t. (5.7) and (5.8) for 
all a G ¥[ti, . . . ,te], and take a K-embcdding t' : ¥[yi, . . . ,yr, zi, . . . , zi] — > ^(IK) 
with r'lf^y^ y^-^ = T s.t. (5.9) for all a G F[ii, . . . , tg]; this is possible by Lemma 5.3. 
Now let s G r(F[2/i,...,y,]) np(F[ii,...,te]), and set / {t')-\s) - t-^{s) G 
¥[yi,...,yr] and g := p~^{s) G F[ti...,te]. To this end, define g' :— (j){g) G 
¥[yi,...,yr,zi,...,zi]. Then by (5.7) we have (5k(5') < (^K(.g)- Since r'(g') = p{g) = 
s and r'(/) = s, and since r' is injective, g' = f. Thus, (5k(/) = S^ig') < S^ig)- O 

Finally, Theorem 5.5 shows that the constructed B G Sum(X) with (5.2) is 
indeed a solution of problem DOS. 

Theorem 5.5. Let (F(yi) . . . {yr),a-) with X = ¥[yi, . . . , j/r] be a UT.^ -field over 
let T : X S{K.) be K-embedding with ev and L, and let (X, ev, (5k) be a sequence 
domain as stated in the General case 4-11; let D := (F(yi) . . . (yr)(si) . . . (su), cr) be 
a Yj^ -extension of {¥[yi) . . . (y,-), (j) where D is a polynomial extension of (F, a) and 
let (5.3) be a R-embedding with a defining function and o-function. Moreover, let 
A G Sum{X) and a G A such that T{a) = (^(fc))fc>o. Then (5K(a) is the Sum{X)- 
optimal depth of A. 

Proof. Take an expression of Sum(X) that produces s := {A{k))k>o from 
a certain point on and that has minimal depth, say d. By Theorem 5.1 we can 
take a E-'-cxtension D (F(?;i) . . . (yr)(ii) ■ • ■ (ie), ct) of (F(?/i) . . . (y^), cr) s.t. D is 
a polynomial extension of (F, cr) and we can assume that there is a K-embedding 
p : X[ti, . . . ,te\ S(K) with a defining function and an o-function with the fol- 
lowing property. There is an a' G . . . ,ie] s.t. p{a') = s and 6-^(0,') < d. By 
Theorem 5.4 (applied twice), (5K(a) — S^{o,')- Moreover, p(a') and r(a) can be 
defined by elements from Sum(X) with depth SKia) by Theorem 4.8. Since d is 
minimal, 5^(0) = d. □ 

6. Application: Simplification of d'Alembertian solutions 

The d'Alembertian solutions (Nor24; AP94; SchOl), a subclass of Liouvillian 
solutions (HS99), of a given recurrence relation are computed by factorizing the re- 
currence into linear right hand factors as much as possible. Given this factorization. 
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one can read of the d'Alembertian solutions which are of the form 

n fc2 fcs-l 

(6.1) h{n) bi{ki) £ 52(^2)- •• ^-(^^) 

for lower bounds ci, . . . , S N; here the bi{ki) and h{n) are given by the objects 
form the coefficients of the recurrence or by products over such elements. Note 
that such solutions can be represented in XlS'^-fields if the occurring products can 
be rephrased accordingly in H-cxtensions. Then applying our refined algorithms to 
such solutions (6.1), we can find sum representations with minimal nested depth. 
Typical examples can be found, e.g., in (DPSW06; Sch07; MS07; BBKS08; OS09; 
BKKS09a; BKKS09b). 

In the following we present two examples with detailed computation steps that 
have been provided by the summation package Sigma. 



6.1. An example from particle physics. In massive higher order calcula- 
tions of Feynman diagrams (BBKS07) the following task of simplification arose. 
Find an alternative sum expression of the definite sum 

(6.2) Sin) 

such that the parameter n does not occur inside of any summation quantifier and 
such that the arising sums are as simple as possible. In order to accomplish this 
task, Sigma computes in a first step the recurrence relation 

- (n + 2){n + if {n^ + 7n+ 16) S{n) 

+ (n + 2) (5n'5 + 62n^ + 318?!.^ + 8Un^ + 1045?i + 540) S{n + 1) 
- 2 (5n^ + Mn^ + 603n* + 2354n^ + 527071^ + 6430n + 3350) S{n + 2) 
+ 2 (5n^ + 96n^ + 78371-* + 34 7877^ + 8906n^ + 12530?! + 7610) S{n + 3) 
- (n + 4) (5n^ + 88n'' + 630?!^ + 2318^^ + 445371 + 3642) S'(77 + 4) 

+ (71 + 4) (71 + 5)'^ (71^ + 571 + 10) 5(71 + 5) = l^ow 

' (n + 3)(n + 4) 

2 (2n^ + 3571^ + 235?i^ + 71871* + 82471^ - 2837i2 - 8697i + lO) 
(71 + l)(7i + 2){n + 3)2(71 + 4)2(71 + 5) 

by a generalized version (Sch07) of Zeilberger's creative telescoping (Zei91). Given 
this recurrence, Sigma computes the d'Alembertian solutions 

Ai = l, A2 = Hn, As^Hl 



^^^(2fc-3)(2fc-l) 

" U ~ i)j 
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^5 = E 



^ 2(fc - 2)(fc - l)fci7fc - (2/s - 1) (3fc2 - 6fc + 2) 
^ ^^^"^^^3 (fc-2)(fc-l)fc(2fc-3)(2fc-l) 



3=3 



ij - 1).? 



i=3 



^ (21-3) {P -3l + 6)B{l) 
'h (2fc-3)(2fc-l) 



where 

I 



~ 2(2r'^-27r'^+117r*-254r^+398r^+2(r-3)(r-2)(r-l)(r+2)_ffrr-446r+204) 
-° ^ {r-2)(r-l)r(r^-5r+10)(r^-3r+6) ' 

r-3 

To be more precise, Ai e Sum„(Q(x)), 1 < i < 5, are the five hnearly independent 
solutions of the homogeneous version of the recurrence, and B G Sum„((!2(x)) is 
one particular solution of the recurrence itself; the depths of Ai , . . . , A5 , B are 
0, 2, 2, 4, 5, 7, respectively. As a consequence, we obtain the general solution 

(6.3) G:=B + ciAi + C2A2 + C3A3 + C4A4 + C5A5 

for constants c^. Checking initial values shows that wc have to choose^° 

ci = ^Cl C2 = j^(48C3 - 67), C3 = §, C4 = i(23 - 8C2), cg = -i 

in order to match (6.3) with S{n) = G{n) for all n e N. 

Finally, Sigma simplifies the derived expressions further and finds sum represen- 
tations with minimal nested depth (see problem DOS). Following the approach de- 
scribed in the previous sections, it computes the IIE'^-field (Q(a:)(/i)(/i2)(ft.4)(i?), cr) 
with a{x) ~ X + 1 and 

a{h) =h + a{h2) = /i2 + j^^, cr(/i4) ^ + j^^, cr(i7) ^ H+ j^^; 

in addition, it delivers a (Q-embedding t : Q{x)[h, h2, ft.4, H] S{Q) with the defin- 
ing function ev : Q{x)[h, /i2, /i4, i/] x N ^ Q canonically given by (2.1) for all / G 
Q{x) and by 

n „ 

ey{h,n)=H^, e^{h2,n) ^ H^^\ ev(/i4, n) = ev(77, n) = ^ 



Moreover, it finds 



ai = 1, a2 = /i, 03 = /i^, a4 = 2 ('''2 ^ '^^) > '^5=2 (^^^ ^ '^'''^ h — 
b^-^{h^ - 48hH + 128/1 - Uhj + {Uh - 69) /12 - 12/14) 



10 



Cfe denotes the Riemann zeta function at k; e.g., 1J2 = 7r^/6. 
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such that ev'{ai,n) = Ai{n) for 1 < i < 5 and ev'{b,n) = B{n). These computa- 
tions lead to the foUowing identities: For ti > 0, 
(6.4) 



A,{n) = ^(^-Hf, + 2H^^^Hn-Hn 



fe=i 



in particular, due to Theorem 5.5, the sum expressions on the right hands sides 
of (6.4) have the minimal depths 0, 2, 2, 2, 2, 3, respectively. To this end, we obtain 
the following identity (BBKS07, equ. 3.14): for n > 0, 



i=i ^ ' A:=l 

6.2. A no ntrivial harmonic sum identity. We look for an indefinite nested 
sum representation of the definite sum 

(6.5) ^(")-E( J 

fc=0 



for similar problems see (DPSW06). First, Sigma finds with creative telescoping 
the recurrence relation 

8(n + l)(2n + 1)^ (64n^ + 480n^ + 1332n^ + 1621n + 735) S{n) - 4(768n* + 8832n^ 

+ 4305671** + 115708n^ + 186452n'' + 183201?i^ + 106442n^ + 33460?i + 4533)S(n + 1) 

+2(n+2)(384n'^+4224n'^ + 18968n^+44610n''+58679n^+42775n^ + 16084n+2616)S(7i+2) 

- (n + 2)(n + 3)^ (64n^ + 224n^ + 276n^ + 141n + 30) S{n + 3) = 

- 3 (576?i® + 4896?!^ + 16660n* + 28761n^ + 26171?i^ + 11574n + 1854) . 

Solving the recurrence in terms of d'Alembertian solutions and checking initial 
values yield the identity 

S{n) = ^ j {-A^{n) - -A,_[n) + B[n)) Vn > 

64j4 - 288 j3 + 468j2 - 323i + 84 



with 



^ " 4z-3 , (j-l).7(2j--3)(4,-7)(4,-3) 



^z(2z-l)' ^ ^[2^-l) 

" ^ ' ^ ^ (j - 1)J(2J - 3)(4j - 7)(4j - 3) 

5 = -E — 



i=2 

j 



i{2i-l) 



where 

j 

\ - 3(2* 

fc(64fc4_544fe3_|_i7i6fc2_2379/c+1227)(64fc*-288fe3+468fc2_323fe+84)(^^'=) 



^ 3(2fc-3)(2fe-l)(4fe-7)(576fc'^-5472fe^+20980fc''-41559fc^+44882fe2-25113fc+5760) 



fc=l 
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So far, this alternative sum representation of (6.5) might not be considered as re- 
ally convincing. For further simplifications, we construct the XlS'^-field {Q{x){b), a) 
with a{x) = X + 1 and a{b) = 2(2x+i) ^' ^^'^ take the Q-monomorphism 
r' : Qix)[b] S{Q) from Example 4.4 (K = Q). Note that the sums Ai,A2,B e 
Sum„(Q(x)[6]) have the depths 2, 3, 5, respectively. Finally, activating our machin- 
ery in this setting (wc represent the sums in a HE'^-field and reinterpret the result 
by an appropriate Q-monomorphism), we arrive at the following identities: For 
n>0, 

A2{n) ^2 (ahI + 4iJ„ + i/|„ + (-4//„ - 2) iJ2„ - i/^l'^) , 

3 " 1 \ 

B{n) =- {aAHI + 16F„ + UHI - (44i?„ + 8) H^n - HH^'J + 14 . 

By Theorem 5.5 we have solved problem DOS for the expressions Ai, A2 and B, 
i.e., we found sum representations with optimal depths 2, 2, 3, respectively. Finally, 
this leads to the following identity: for n > 0, 

(';) = (';) (4^- - 4i/2„i?„ + - Hi^j + ' ^ ^) ■ 
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